Local physics of magnetization plateaux in the Shastry-Sutherland model 
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We address the physical mechanism responsible for the emergence of magnetization plateaux in 
the Shastry-Sutherland model. By using a hierarchical mean-field approach we demonstrate that a 
plateau is stabilized in a certain spin pattern, satisfying local commensurability conditions derived 
from our formalism. Our results provide evidence in favor of a robust local physics nature of the 
plateaux states, and are in agreement with recent NMR experiments on SrCu2(B03)2. 
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Introduction.- The interplay between quantum me- 
chanics and the atomic lattice topology often leads 
to a complex mosaic of physical phenomena in low- 
dimensional frustrated magnets [ij. A prominent rep- 
resentative of this class of materials is the layered com- 
pound SrCu2 (603)2, "^hich recently received a lot of at- 
tention because of its fascinating properties in an exter- 
nal magnetic field h, namely the emergence of magnetic 
plateaux at certain fractions of the saturated magneti- 
zation Msat • The first experimental observations of the 
plateaux were reported in Q for m = M/Mgat = 1/8 and 
1/4, and somewhat later for m = 1/3 Subsequent nu- 
clear magnetic resonance (NMR) experiments |3, IH re- 
vealed spontaneous breaking of the lattice translational 
symmetry within the 1/8 plateau, and also indicated that 
the spin superlattice persists right above this fraction 6 1 . 
The field was reignited by the work of Sebastian et al. [7| , 
where additional plateaux at exotic values m = 1/9, 1/7, 
1/5 and 2/9 were reported. However, direct observation 
of the emerging spin superstructures remains an exper- 
imental challenge, primarily due to the high magnetic 
fields (~ 30 — 50 Tesla) involved in measurements. 

The nature of the magnetic states and physical mech- 
anism leading to the plateaux are also yet to be under- 
stood. It is believed that the Heisenberg antiferromag- 
netic model on a frustrated Shastry-Sutherland (SS) lat- 
tice with N sites [i| (Fig. [J), 




h = jY,s,-s, 



2JaY,S,-Sj-hY,S!, 



(1) 



captures the essential magnetic properties of 
SrCu2 (603)2 relatively high fields. In Eq. (P) 
Si denotes a spin-1/2 operator at site i; the first sum is 
the usual nearest-neighbor (NN) Heisenberg term, while 
the second one runs over dimers; J and a ^ 0. This 
model is quasi-exactly solvable @ for a ^ 1 + h/2: the 
ground state (GS) is a direct product of singlet dimer 
states, and was shown to be stable up to a 0.71-0.75 
in zero field In general, it is an intractable quantum 
many-body problem where approximation schemes are 
needed to deal with large- A'' systems. 

All theories proposed to address this unusual magne- 



FIG. 1: The SS lattice. Circles denote spins, dashed lines 
correspond to the NN J coupling, double solid lines denote 
next-NN interactions (dimers). The simplest choice of a de- 
gree of freedom, which does not cut dimers, is shown in black. 



tization phenomenon start from the SS model. However, 
the physical mechanism stabilizing the plateau states, 
their nature, and the structure of the magnetization 
curve are still actively debated. Current ideas can be 
broadly divided into two groups. The first one advocates 
subtle non-local (in the spins) correlations leading to an 
underlying spin structure which preserves lattice symme- 
tries [lO[ . Technically, it employs a mapping of the orig- 
inal spins to fermions coupled to a Chern-Simons gauge 
field, and then performs a Hartree-Fock decoupling. In 
this way, the qualitative shape of the SrCu2 (603)2 ™^^S" 
netization curve was reproduced in high fields, but the 
lowest plateau at 1/8 was missing. Later, this non-local 
mean-field approach was extended to include inhomoge- 
neous phases 0, and it was argued that the plateaux 
correspond to stripe states with broken lattice symme- 
tries. Remarkably, the length scale ^ associated with the 
emerging spin superlattice was found to be ^ ~ 100 lat- 
tice spacings. The second group contends that the mag- 
netization process can be described in terms of polarized 
dimers (triplons), which propagate in the background of 
singlet dimers [3, [HI. They developed effective hard- 
core boson models (truncating the original dimer Hilbert 
space), solved by perturbative [l3| or CORE [l3| tech- 
niques, and found that the plateaux states correspond to 
crystal phases with ^ ~ 10 lattice constants. 
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Such diversity of theoretical predictions demands fur- 
ther investigation. In this paper we use the hierarchical 
mean- field (HMF) method IJ], in an attempt to clarify 
the nature and physical mechanism, responsible for the 
emergence of magnetic plateaux in the SS model. Un- 
like previous calculations, we deal directly with the SS 
Hamiltonian ([T]), not with effective Hamiltonians as in 
Refs. [H, and combine exact diagonalization data 
with a simple and controlled approximation for the GS 
wavefunction. For instance, we do not discard the M = 
dimcr triplet states, necessary for the propagation of a 
triplon. We focus on higher-lying fractions, whose ex- 
istence has been confirmed experimentally. Our results 
support the local physics nature of the plateau states. In 
particular, it is explicitly demonstrated how to construct 
those robust states based on a set of commensurability 
rules that we derived. Our conclusions are also in agree- 
ment with the interpretation of NMR measurements Q . 

Method. - The HMF approach is based on the assump- 
tion that the physics of the problem is local in a partic- 
ular representation. Since the SS model is formulated 
in terms of localized spins, it is natural to work in real 
space. The main idea of our method revolves around the 
concept of a relevant degree of freedom - a spin cluster - 
which is used to build up the system. Essential quantum 
correlations, which drive the physics of the problem, are 
captured by this local representation. The SS Hamilto- 
nian is then rewritten in terms of these coarse-grained 
variables and a mean-field decoupling is eventually ap- 
plied to determine properties of the system. The method 
is only limited by finite-size effects and becomes asymp- 
totically exact in the thermodynamic limit. Thus, the 
(generally) exponentially hard problem of determining 
the GS of the model is reduced to a polynomially com- 
plex one. By identifying each state of a cluster with a 
Schwinger boson (SB) and computing matrix elements 
of spin operators between these states, one can rewrite 
exactly the Hamiltonian of Eq. ([T]) in terms of new (clus- 
ter) variables 



H 



Here the repeated indices a, b, etc., which label states 
of an iVq-spin cluster, are summed over, and i denotes 
sites in the coarse-grained lattice. The operators jj^ 
that create a particular state of a cluster are SU(2^''^ 

SBs subject to the constraint jjalia = 1 on each 
site; ea sltc exact cluster eigenenergics. Since the orig- 
inal SS Hamiltonian involves only two-spin interactions, 
in the new representation there will be only two-boson 
scattering processes: the corresponding matrix elements 

a' h' 

are denoted by [Hf^^)^^ . The second term in Eq. ([2|) 
describes the renormalization of the cluster energy due 
to its interaction with environment. Thus, our method 
deals with an infinite system and finite-size effects enter 



only through a particular choice of the cluster. The sym- 
bol (ij)cr, with cr = 1, 2, . . ., indicates pairs of neighboring 
blocks, coupled by the same number of J-links. 

Application of the HMF method to the SS model starts 
by recalling that the phases within plateaux break the 
lattice translational invariance. Therefore, the best solu- 
tion will be obtained, if the degree of freedom matches 
the unit cell of the spin superstructure. For each cluster 
size, Nq, and magnetization 
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we determine the lowest-energy configuration (i.e., the 
cluster shape and corresponding tiling of the lattice). By 
virtue of previous argument, this solution will have the 
"right" symmetry. Performing this operation for suc- 
cessive values of Ng up to the largest one that can be 
handled; we obtain a set of magnetization plateaux to- 
gether with their corresponding spin profiles. It follows 
that the particular choice of coarse graining is critical for 
the success of this program. One should recall that the 
experimental value for a is 0.74-0.84, i.e. the intradimer 
coupling seems to be "more relevant" than the interdimer 
one. Therefore, it is natural to consider only those clus- 
ters, which contain an integer number of a-links. This 
constraint turns out to be quite severe. It follows that the 
degree of freedom must also contain an integer number 
of "minimal" blocks, shown in Fig. [T] in black: other- 
wise the tiling of the lattice will not be complete. These 
requirements comprise a set of local commensurability 
conditions, necessary to stabilize a plateau. 

Another crucial issue is the way the interaction terms 
in Eq. ([2]) are handled. In an attempt to simplify mat- 
ters, we use the straightforward Hartree approximation, 
i.e., we consider the trial GS wavefunction 



(3) 



Here |0} is the SB vacuum and Ra are variational param- 
eters, which constitute the cluster wavefunction. Since H 
is real-valued, we can choose Ra to be also real. Clearly, 
this state is cluster translationally invariant and has ex- 
actly one boson per coarse-grained lattice site (so the 
constraint is exactly satisfied). Next, we compute the 
expectation value of H in the state ([3]), subject to peri- 
odic boundary conditions, and minimize it with respect 
to Ra- In this manner one obtains the approximate GS 
energy Eq as a function of the magnetic field h. 

It is important to emphasize the simplicity of our ap- 
proach. By using a more sophisticated ansatz (e.g. a 
Jastrow-type correlated wavefunction 14| ) , we could im- 
prove energies but the physical mechanism and robust 
structure of the plateaux will remain intact. Despite its 
simplicity, the ansatz of Eq. ^ was accurate enough to 
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TABLE I: Representative values of the GS energy parameter 
eo- Numbers in parenthesis denote the size of a cluster, Nq. 
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FIG. 2: Schematic spin profiles within plateaux. (Gray) black 
circles correspond to polarizations (anti) parallel to the field; 
their sizes encode the magnitude of the local magnetic mo- 
ment. Empty circles denote sites with \{Si)\ < 10~^. The 
clusters used in HMF calculations consist of dark gray dimers. 
Light gray dimers represent the NN cluster. For m = 1/3, 1/4 
and 1/6 dark and hatched dimers constitute the 24-spin clus- 
ter. Thin lines indicate unit cells of the spin superlattice. 



yield the quantitatively correct phase diagram of the Ji- 
J2 model which involves gapless phases. In the SS 
model states within the plateaux are gapped, therefore, 
our method should be tailored for this problem. 

Results. - To guarantee that our results reduce to the 
exact solution in the limit h 0, we mainly consider 
the region a > 1. The simplest degree of freedom, con- 
sisting of 4 spins, is shown in Fig. [TJ Using this cluster 
in our HMF scheme one obtains stable plateaux only at 
TO = 1/2 and m = 1. Clearly, larger blocks are neces- 
sary to stabilize plateaux at lower magnetization frac- 
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tions. 

fc = 2, 



Here, we consider cluster sizes iV„ 



4fc with 



,6 and discuss only plateaux at 1/3, 1/4, 1/5, 
1/6 and 1/8, supported in minimal clusters of Nq = 12, 



8, 20, 12 and 16 spins, respectively {Mgat = Nq/2). In 
Fig. [2]we present local spin profiles, corresponding to the 
lowest energy configurations, for each of these fractions. 
Comparison of patterns for different plateaux, shows that 
states 1/n with n even are characterized by one polar- 
ized dimer per unit cell, while cells of odd-n states have 
two triplons. For a given plateau, there typically exist 
several possible coarse-graining scenarios, characterized 
by different clusters and tessellations of the SS lattice, 
but identical unit cells. Although these configurations 
have slightly different energies, their existence provides 
an important check for robustness of local correlations, 
stabilizing the plateau states. The patterns for all frac- 
tions except 1/5, are similar to those obtained in Refs. 
[1, Strictly speaking, the profiles in Fig. [2] are well 
defined only for large values of a > 1 and quickly smear 
out with decreasing a. This effect is difficult to capture 
within the effective model calculations, like [l3 |. 

A clear advantage of our approach, compared to pre- 
vious works, is its ability to compute GS energies of 
the original SS model. Within each plateau we have: 
EQ{h)/N = — £0 — mh/2. The parameter eg is presented 
in Tablc[T]for some values of a and different cluster sizes. 
In order to address finite-size effects, in Fig. |3]we present 
the high magnetic field phase diagram of the SS model for 
a ^ 1. All fractions were calculated using the largest pos- 
sible cluster. Due to the insulating nature of the plateau 
states, the finite-size corrections are not expected to sig- 
nificantly affect their stability. For instance, for to = 1/6 
the energy difference between 12- and 24-spin clusters is 
only 5% of its width for the values of a shown in Ta- 
ble m This observation serves as additional evidence in 
favor of a universal physical mechanism leading to the 
plateaux. 

As it was already mentioned, the HMF method does 
not involve truncation of the dimer Hilbert space. In 
order to understand consequences of this approximation, 
we computed Eq for different plateaux, ignoring the dimer 
state III). The resulting absolute error is of the same or- 
der of magnitude as finite size effects and plateau widths 
(cf. Table [H, which leads to a sizable change in the rela- 
tive stability of the plateaux. For example, at a = 1.1 the 
average error is 3 • 10"'^ J, and the lower boundary of the 
1/8 state shifts by 0.04 J. Therefore, conclusions of the 
effective boson model calculations, which employ similar 
truncation, should generally be taken with caution. 



4 



o.ob 

0.08 
0.07 
"5 0.06 
g 0.05 
4" 0.04 
-£- 0.03 
0.02 
0.01 






1/4 

1/5 // 


- 


/ 1/6 







0.2 0.4 0.6 0.8 


1/3" 






1/4 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

1/a 

FIG. 3: High magnetic field phase diagram of the SS model 
for a ^ 1. ho{a) denotes the field after which the first plateau 
(at 1/8) emerges. Fractions indicate values of m. For q 3> 1 
the triplons (•=•) in Fig. [2] become fully polarized and other 
dimers within the clusters turn into perfect singlets. 

Discussion.- Although the effective model approach 
does yield a sequence of plateaux, their understanding 
remains incomplete. Our work addresses this issue by fo- 
cusing on the nature and correlations of the magnetic 
plateau states. In particular, the analysis presented 
above, allows the formulation of a set of universal rules 
leading to well-defined spin patterns (Fig. which can 
be probed, e.g. by neutrons. These rules define a hierar- 
chy of variational plateau wavefunctions and constitute 
a central prediction of our work. For a robust state to 
emerge at a given magnetization fraction ?7i, the commen- 
surability conditions that have to be fulfilled are: (i) the 
(cluster) degree of freedom must contain an even number 
of dimers; {ii) the SS lattice must be tessellated com- 
pletely with these clusters; {iii) the size of the cluster 
(unit cell), Nq, must allow the plateau state at m, there- 
fore, Nq = 2M/m with M = l,...,Nq/2 chosen in a 
way such that Nq is divisible by four; (iv) the number 
of triplons (^— ^) per cluster is M and its shape must 
be such that each triplon is surrounded by two dimers of 
the type =•, within this cluster. The application of the 
above constraints leaves us with the essentially combina- 
torial problem of actually determining the symmetry and 
periodicity of the spin superstructure (see Fig. [J). 

There also exists a number of concrete discrepancies 
between our work and recent publications [l^ . [isj , which, 
nevertheless, support our general conclusion regarding 
the local nature of the plateau states. First, the exper- 
imentally observed plateaux at 1/4 and 1/8, which we 
found to be quite robust, were claimed to be unstable 
in However, the magnetization profile, presented in 
Fig. [2]for m — 1/8, which persists at a = 0.787, adequate 
for SrCu2(B03)2, is consistent with the interpretation of 



available NMR data [1, Q for this material. We believe 
that the origin of these states is purely magnetic and no 
additional interactions beyond the SS model are required, 
in contradiction with the claim of Ref. [l3l ■ Our results 
also yield a stable 1/5 plateau, contrary to the conclu- 
sions of Refs. [ll, [l3| . We note that this fraction was ob- 
served in torque measurements of 0, however, their pro- 
posed spin superlattice differs dramatically from the one 
predicted in our Fig. [21 Another distinction concerns the 
robustness of the 1/6 plateau advocated in [l^l, which, 
although present in our calculation, has a significantly 
smaller relative width (see the discussion above). Other 
fractions at 1/9, 2/9 and 2/15, observed in Refs. and 
[l3| . can also be obtained within our approach, but this 
requires significantly larger clusters than the ones used 
here. By virtue of our commensurability arguments, we 
expect the plateaux at 1/9 and 2/9 to emerge in degrees 
of freedom containing at least 36 spins, while the 2/15 
fraction will be stabilized in a 60-spin cluster. 

Finally we note that the precise shape of the magne- 
tization curve (the relative energy stability of different 
plateaux) is quite sensitive to the value of a (i.e., the 
particular compound) and, most importantly, since there 
is no exact solution of the SS model at these high fields, 
it depends on the particular approximation scheme. Ex- 
perimentally, other physical interactions not included in 
the SS model may also add to this uncertainty. 
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